Observational constraints on a unified dark matter and dark energy model 

based on generalized Chaplygin gas 
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We study a generalized version of Chaplygin gas as unified model of dark matter and dark energy. 
Using realistic theoretical models and the currently available observational data from the age of 
the universe, the expansion history based on the type la supernovae, the matter power spectrum, 
the cosmic microwave background radiation anisotropy power spectra, and the perturbation growth 
factor we put the unified model under observational test. As the model has only two free parameters 
in the flat Friedmann background [ACDM (cold dark matter) model has only one free parameter] 
we show that the model is already tightly constrained by currently available observations. The only 
parameter space extremely close to the ACDM model is allowed in this unified model. 

PACS numbers: 98.80.-k, 95.36.+X 
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I. INTRODUCTION 



A unified dark matter and dark energy model moti- 
vated by the (generalized) Chaplygin gas has been intro- 
duced in the literature P, d, Q , and has been extensive^ 
investigated in the cosmological studies 
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gated m the cosmolo gica l studies I 




Many constraints on the Chaplygin gas model parame- 
ters have been placed (predicted) based on current (up- 
coming) astronomical observations such as typ e la super- 
novae (SNiaWii, pj3Ji, [il p, [H, 113, H H M, m, 

f.j2a. [26l l27l. [28l. [29l. Isil. l33l. [35| . large-scale structure 
0, M, ^M, [23, [Sllj [33j [38 J4p| , cosmic microwave back- 
ground radiation jX^MBTjini i, 0, IM S , gravi- 
tational Icnsing [Til. IT^. [Tsl [20| . gamma-ray bursts [sot . 
X-ray luminosity of galaxy clusters [T^, ,22] , look-back 
time-redshift data [Hjlsl], angular size-redshift data 24 1, 
Hubble parameter-redshift data [13, UHl , FanarofF-Riley 
type lib radio galaxies [l3;[ll|, and so on. 



It is known that such a unified model is constrained 
by observations so that parameter space close to the con- 
ventional ACDM models is allowed. But current status 
of the model in parameter space different from close-to- 
ACDM has been unclear. The distancc-redshift relation, 
the age of the universe, abundances of light elements, 
the large-scale matter density and velocity power spec- 
trum, the CMB temperature and polarization anisotropy 
power spectra, and the perturbation growth factor can be 
regarded as the main pillars of modern cosmology where 
theories meet with observations. Here, we investigate the 
unified model by comparing realistic theoretical predic- 
tions with currently available observations. Our conclu- 
sions can be found in section HVl 



II. GENERALIZED CHAPLYGIN GAS 



We introduce a fluid X with an equation of state 



Px 



(1) 



where p and // are the pressure and the energy density, 
and A and a are constants; tildes indicate covariant quan- 
tities. The Chaplygin gas is a case with a = 1. It is in- 
teresting to note that any fluid with barotropic equation 
of state, px = pxifi-x), has an exact scalar field theory 
counterpart with a tachyonic kinetic term given by an 
action (we set SttG = 1 = c) 



S 



X (X e-' i^x+px , (2) 



where X = ^■'^(p^c, see Appendix A for a proof. For a 
fiuid in Eq. (QJ we have 



Px 



1-[XIB 



(3) 



where S is a constant. We call the fluid or field based on 
Eqs. ([T]) and ([3]), a generalized Chaplygin gas (GCG). 

To the background order in Friedmann world model 
the energy conservation equation for GCG component 
gives 



Hx 



wx 



A 



l+a _ 1 + 



a3(l+a) 



px 



px = 
-1 



a3(i+") 



2 px 

Cx = — = -awx, 

MX 



(4) 



where a is the scale factor and A*^" > A (we set 
ao = 1 at the present epoch). The equation of state 
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parameter wx is negligible in the early era thus the 
GCG potentially acts as a dark matter, and approaches 
wxo ~ ^^/Ma'o"(-^ ^1) recent past thus it poten- 
tially acts as a dark energy depending on suitable choice 
of parameters. 

We consider a flat background with baryon, radiation 
(photons and neutrinos), and the GCG. At the present 
epoch the Friedmann equation gives 



xo 



1 - n 



bO 



(5) 



where fli = fii/ fJ-crit with ficrit the critical density. Thus, 
/ixo is completely fixed by the Friedmann equation. 
Equation Q gives 



(6) 



As A is determined by wxo and a, in the following we 
will regard wxo and a as the two free parameters to be 
constrained by comparing the theoretical consequences 
with observations. 

The ACDM limit can be reached by taking a = 0, and 
identifying A = = A and fico = Mjco — A; indicates 
density of the CDM. Thus, in a flat background, from 
Eqs. ([5|) and ([6]) we have 



WXG = 



^AO 

n 



XO 



cO 



(7) 



We can show that even linear perturbation system (i.e., 
equations for dc and Vc) also coincides exactly with the 
ACDM case. As a fiducial model, we use flat ACDM 
model consistent with Wilkinson Microwave Anisotropy 
Probe (WMAP) 5-year data [4l| {^to = 0.0456, ^co = 
0.2284, Oao = 0.726, h = 0.705, n, = 0.960, as = 0.812, 
To = 2.725 K, Fho = 0.24, and = 3.04), but 
without reionization history, see [i^. Thus we have 
Wxo = -0.7607. 

In Fig. [T] we present the likelihood contours in a-wxo 
plane based on the SNIa data. Also presented are the 
age lines. Notice that the likelihood contour as well as 
the age lines depend on flbo in our GCG model. It is 
interesting to see that, except for small region close to 
the ACDM, most of the la region based on SNIa can 
be ruled out by applying age of the universe larger than 
13.5 Gyrs. Notice that the SNIa alone favors a = 0.9971 
and Wxo = —0.8523 (indicated by a bullet) which will be 
completely excluded by the matter and the CMB power 
spectra (see Fig. [3]) and also by the age limit. Background 
evolutions for the locations indicated by * in this Figure 
are presented in Fig. [2] Perturbation power spectra and 
the growth factors for the same parameters are presented 
in Figs. [3] and H 

We consider scalar-type perturbations. For a fluid with 
barotropic equation of state we naturally have 



Spx = CxSfix, 



(8) 
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FIG. 1: Likelihood (l-3cr) contours in a-wxo plane based 
on the SNIa Union sample The locations of ACDM 

(a = 0) based on the WMAP (Oao = 0.726) and Union 
data (r^Ao ~ 0.713) are indicated as x . The maximum like- 
lihood point is indicated by a bullet (•). Locations indi- 
cated by ★ are studied in details in Figs. [2]|4l Locations near 
Of = at Wxo = —0.7607 are studied in details in Figs. I5I7I 
We marginalized over the Hubble constant. Short-and-Iong- 
dashed lines are age lines in the unit of (/i/0.705)~^ Gyrs. 
The age lines together with the SNIa data already show that 
parameters with q > are largely excluded by currently 
available observations. A short-dashed line shows parame- 
ters which reach c^g = 1 at the present epoch; the right side 
of the line has c^g > 1 thus becoming super-luminal. 



stress also vanishes the GCG is exactly an ideal fluid. 
Although the fluid deflnition of the GCG model is am- 
biguous about the ideal fluid nature, this can be proved 
by the fleld theoretic counterpart based on Eq. ([2]); see 
Appendix A for the proof. Without taking the temporal 
gauge condition the energy and the momentum conser- 
vation equations for the GCG become 



Ox 



Vx 



3 {wx 



^x 



HS 



X 



-|- (1 + Wx) n — 3Ha Vx 

' a 



(9) 



k 



1 - ic'x) Hvx + - 



-X 



<5a + -«,(10) 
a \ + Wx a 



where Cy = —awx for GCG; these are Eqs. (A7) and 
(A8) in |44| with 5x = SfJ.x/fJ'X, vx a perturbed velocity 
of GCG component, a a perturbed metric variable and 
K a perturbed part of the expansion scalar of the normal 
frame, see [47l |: k is the comoving wave number. 



III. OBSERVATIONAL CONSTRAINTS 



without the entropic perturbation. This relation is valid 
without taking any gauge condition. As the anisotropic 



In order to calculate the matter and CMB power spec- 
tra, and the baryon density perturbation growth factor 
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FIG. 2: Top panels: Evolution of il^ and \x.i as a func- 
tion of scale factor ait) for values of (a,wxo) indicated as 
★ in Fig. [T] i = r, b, indicates radiation, baryon, and 
GCG, respectively. Middle and bottom panels: Evolution 
of fix, wjf, Hx{z) / Hi^coM {Hx is the Hubble parame- 
ter in our GCG model), and the relative distance modulus 
Aij,{z) = Hx (z) — /XACDM (z) for the same set of GCG param- 
eters. In all panels, ACDM predictions are shown as thick 
black curves. In the A/x-plot, the grey open squares with er- 
ror bars represent the deviation of SNIa data points from the 
fiducial ACDM model considered here. The binned SNIa data 
are based on the Union sample [4^ . 



we solve a system composed of baryon, radiation (han- 
dled using the Boltzmann equation or tight coupling ap- 
proximation), together with the GCG described by Eqs. 
© and pop representing the dark matter and the dark 
energy in a unified way. We consider a flat background 
with similar parameters as our fiducial ACDM model 
mentioned below Eq. ([7]) including neutrino components. 
Our set o f eq uations and the numerical methods are pre- 
sented in |45| , 

We solved the system in three different gauge con- 
ditions: the synchronous gauge (SG), the uniform- 
expansion gauge (UEG), and the uniform-curvature 
gauge (UCG); see [i^, 113, Ell for the description of the 
gauges. The matter power spectrum and the perturba- 
tion growth factor are for the baryonic matter perturba- 
tion in the synchronous gauge which is a gauge-invariant 
concept in our situation. The CMB temperature and po- 
larization anisotropies are naturally gauge invariant. The 
final results of these gauge-invariant variables calculated 
in our three different gauge conditions should coincide; 
this provides a numerical check of the calculations. 

In Figs. [2] and [3] we present the background evolu- 




FIG. 3: The matter (baryon) power spectrum (top-left), and 
CMB TT (top-right), EE (bottom- left), TE (bottom-right) 
power spectra of GCG models with parameters used in Fig. 
[2l and the same colored code. All calculations are made in 
three different gauge conditions (SG, UEG, and UCG), and 
the results in the three gauges coincide exactly. The mat- 
ter and CMB power spectra of the ACDM model have been 
normalized with as and COBE spectrum, respectively. For 
comparison, all the GCG power spectra have been normal- 
ized with the ACDM ones at fc = 0.01 h Mpc"^ for matter 
power spectrum and £ = 700 for CMB ones. The ratios of 
GCG powers to ACDM predictions are also shown in the bot- 
tom region of top panels. For matter and CMB TT power 
spectra, recent measurements from SDSS DR7 Luminous Red 
Galaxies (LRG) [H and WMAP 5-year [4l[ data (including 
the cosmic variance) have been added (grey dots with error 
bars) together with fractional errors of observed spectra. For 
a correct comparison with the LRG band powers, the model 
power spectrum should include the convolution effect caused 
by LRG band power window functions and the non-linear 
clustering information. 



tion, and the matter and the CMB power spectra for 
several typical parameters indicated as * in Fig. [TJ The 
bottom-right panel in Fig. [5] shows that all the parame- 
ters we consider fit well with the SNIa data. However, 
the power spectra in Fig. [3] show that all the GCG mod- 
els we consider fail to fit the observed matter and CMB 
power spectra simultaneously. The baryon density per- 
turbation growth factor presented in Fig. [4] also confirms 
this result. 

Thus, we find that the only region in a-wxo parameter 
space allowed by the current observations is the location 
close to ACDM. In order to constrain the observationally 
allowed variation of a in that region, in Figs. [5] and [6] 
we investigate the case of — 10~^ < a < 10~^ for fixed 
wxo = —0.7607. In Fig. [5] we present the matter and the 
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FIG. 4: Evolution of baryon density perturbation (top-left), 
and the normalized perturbation growth factor g = ja) in 
three different scales for the same parameters used in Figs. [2] 
andO We add 1% error bar expected from future X-ray and 
weak lensing observations [49| . 
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FIG. 5: The same as in Fig.[3]for values near q = (indicated 
in the Figure) with uixo = —0.7607. 



CMB power spectra; we do not present the background 
evolutions which are quite similar to the ACDM ones in 
Fig. [21 For the parameters considered the CMB power 
spectra are similar to the ACDM, thus observationally in- 
distinguishable for a < 10~^. The matter power spectra 
in Fig. [SI however, depend more sensitively on the value 
of a. For example, for a > 10~^ the current observation 
can be used to distinguish its deviation. The baryon den- 
sity perturbation growth factor presented in Fig. [5| also 
confirms this result which shows that the deviations are 



FIG. 6: The same as in Fig. |3] for the same parameters used 
in Fig. [5] 



TABLE I: GCG model parameter constraints (68.3% CL) 
from SDSS DR7 LRG and ACDM-motivated mock power 
spectrum data based on likelihood distribution around a — 
(close to ACDM model) . 
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particularly significant in the small scale. 

Based on severe oscillations and divergent behaviors in 
the small scale of GCG power spectrum for nonvanishing 
a, authors of [l^l concluded that \a\ larger than 10~^ 
arc excluded by the observation. This is understandable 
because the sound velocity squared becomes nega- 
tive/positive, thus causing instability/oscillation for neg- 
ative/positive a in the small-scale limit, see Appendix 
B for the analysis. However, later it has been shown 
that, despite the heavy oscillation and divergence of the 
GCG power spectrum, the accornpanied baryon power 
spectrum behaves relatively well d, Ull . In order to re- 
solve the issue clearly, in Fig. [7| we present the baryon 
power spectra together with the GCG power spectra for 
the same parameters used in Fig. [S] The Figure con- 
firms that despite the wild oscillations and divergences 
of GCG power spectra, the baryon power spectra behave 
much mildly. The Figure shows, however, that although 
\a\ ~ 10^^ arc surely acceptable, a ~ 10~^ gives devi- 
ation compared with current observation and a ^ 10~^ 
could be already excluded. For negative value of a wc 
have more stringent limit so that a < —5.0 x 10~^ al- 
ready shows diverging behavior in the small scale limit. 

In Figs. [5][7| we show that near a = (thus for GCG 
models close to ACDM) the observed matter power spec- 
trum can be used to constrain tightly the allowed GCG 
models, whereas the CMB power spectra are relatively 
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FIG. 7; The baryonic matter power spectrum in Fig. (5] 
together with the GCG power spectrum plotted in dashed 
lines. Although the baryon power spectra are normalized at 
k = 0.01 h Mpc~^, the GCG power spectra are normalized at 
k = 0.001 h Mpc~^ due to heavy oscillation of the latter ones. 
Compared with the GCG power spectra which oscillate signif- 
icantly, baryon power spectra relatively behave mildly. The 
grey and blue dots with error bars represent the SDSS DR7 
LRG and the ACDM mock power spectrum data, respectively. 
For the latter, we define 50 data points in the logarithmic in- 
terval between 0.02 and 0.3 /iMpc~^, and perturb the linear 
ACDM matter power spectrum at each point by adding a 
Gaussian noise with 10% of the power spectrum amplitude as 
standard deviation. This ACDM-motivated mock data will 
be used to constrain our GCG model in Fig. |5]and Table I. 



SHTiilar to ACDM. In Fig. [5] we show the hkelihood con- 
tours of GCG model parameters obtained from the SDSS 
DR7 LRG and ACDM-motivated mock power spectrum 
data in Fig. [7] (see Appendix C for the method). Table U 
lists 68.3% confidence limits (CL) of GCG model param- 
eters estimated based on likelihood distribution around 
the narrow region close to ACDM model (inner panel of 
Fig. [5]). Although not favored in the age constraint, we 
find another (besides the close-to-ACDM) region in the 
ct-wxo parameter plane where the matter power spectra 
are favored by current observations. In Fig.[9]we show the 
matter and the CMB power spectra for several parame- 
ters within la domain in the matter power spectrum in 
Fig. [5] Authors of |38| showed that for large a with wxo 
extremely close to —1 the matter power spectra are com- 
patible with observation; this is consistent with our Fig. 
[5] (note that the truncation of contour at a ~ 4.6 is due 
to the finite bin size of wxo used in the GCG model pa- 
rameter search). However, our Fig. [9] shows that despite 
the observational success in the matter power spectra, 
significant deviations in the CMB power spectra are un- 
avoidable, thus this additional domain is excluded by the 
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FIG. 8: Likelihood (l-3a) contours in a-wxo plane based on 
the matter power spectrum measurements from SDSS DR7 
LRG [4^ (red, yellow, green curves). For the method used, 
see Appendix C. Notice that besides the region around ACDM 
near a = (see the inner panel), the matter power spectrum 
favors another island with positive a. This island, however, 
can be excluded by the CMB power spectra in Fig. [9] and 
the age limit in Fig. [1] In the both panels, we add similar 
constraint based on ACDM-motivated mock power spectrum 
data in Fig. [7] (blue contours). GCG model parameters used 
in Figs. [2]|4] are indicated by power spectra for several pa- 
rameters consistent with LRG power spectrum indicated by 
+ are presented in Fig. |9l 
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FIG. 9: The same as in Fig. |3] for parameters indicated. 
We present parameters which show similar behavior as the 
ACDM in the matter power spectra. All these models, how- 
ever, show severe deviations in the CMB power spectra, thus 
are excluded. 
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CMB observations (in addition to the age test). 



IV. DISCUSSION 

We have systematically studied observational conse- 
quences of GCG unified dark matter and dark energy 
model. In the background Friedmann world model the 
basic requirement of minimum age of the universe to- 
gether with the SNIa luminosity-redshift data already 
leaves only small parameter space close to the ACDM 
model, see Fig. [TJ In the perturbation study, we have 
shown that the matter power spectrum and the CMB 
power spectra are mutually exclusive except for narrow 
region close to the ACDM model, see Figs. [31 [51 and [HI 
A tight constraint on a-paramcter in that allowed re- 
gion can be obtained from the (baryonic) matter power 
spectrum as --5 x 10^^ < a < 10^'*, see Figs. [5l and 
[71 More realistic constraint can be found in Fig. [HI and 
Table [J as |a| < 10'^ from the SDSS DR7 LRG power 
spectrum and \a\ < 10~^ from the ACDM mock power 
spectrum data. Although this allowed region is wider 
than what was concluded based on the GCG power spec- 
trum in (see Fig. \7}, the constraint based on the 
baryon power spectrum is still severe enough so that the 
observationally allowed GCG model can be regarded as 
extremely close to the ACDM model, thus, effectively in- 
distinguishable from the ACDM model. Although named 
as a unified model the GCG model has two free param- 
eters, a and wxo, whereas the ACDM has only one free 
parameter 17 aq. From the perspective of more wider pa- 
rameter space of our GCG model, it is remarkable to 
notice the distinguished success and its uniqueness of the 
ACDM model. 
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Appendices 

Appendix A. Tachyonic field correspondence: 

Action in Eq. gives 

Tab = Pxgab - PX,x4>,a4>,b- (H) 

Under the energy frame, with qa = 0, we have 

11= -pX +'2.pX.xX, P = PX, TTah = 0, (12) 



where the fluid quantities are defined in Eq. (2) of [5C 
From px = —px + '2px,xX we can derive 



(13) 



Thus, for any barotropic fluid with px = px{fix) we 
have a corresponding tachyonic field with px{X) given 
by the above relation. To the background order, we have 
c\ = pxl px = Px,x I lix,x- To the perturbed order, we 
have 5px = Px,x^X and biix = ^ix,x^X, thus ex = 
5px — c\5^x = 0. As we have ex = and irxab = 0, the 
tachyonic field based on Eq. ([2]) corresponds to an ideal 
fluid. 

Appendix B. Small-scale instability for ct < 0: 

For a < we have < 0. As the GCG is an ideal 
fluid cx can be interpreted as the sound velocity, and 
imaginary sound speed naturally leads to small scale in- 
stability. In order to show the instability we take the 
GCG-comoving gauge which sets vx = 0. Together with 
the Raychaudhury equation (see Eq. (A14) in [i^ D 

k + 2Hn = - ( 3i/ - a + 47rG {Sfi + Sdp) , (14) 



Eqs. dH) and ^ lead to 

Sx + {2- 6wx 4 ''"^ 
- 3 (5wx - 3cx) 



34) HSx 



3 {wx 



-X. 



H 



^x 



1 + Wx 

+At:G (1 



X +Pxj 



H 



w 



[l + 3c;) p5x 



4ttG (1 + wx) {S^l + 3Sp) ~ cjc—Sx, 



(15) 



where Sfj, and Sp are collective (total) perturbed energy 
density and perturbed pressure; w = p/p, and = p/fi 
where ^ and p are collective (total) energy density and 
pressure of the background world model. We note that 
the above equation is generally valid in the presence of 
other components (baryon, radiation, etc.) and the back- 
groimd curvature. The presence of Cxk^Sx term in the 
right-hand-side of above equation shows that for nega- 
tive/positive we have strong pressure-caused instabil- 
ity/oscillation in the small scale limit where k is large. 

Appendix C. Likelihood estimation of GCG 
model parameters using the matter power spec- 
trum data: Here we briefly summarize how the like- 
lihood contours of GCG model parameters in Fig. [51 
have been obtained from the observed and mock power 
spectrum data. The power spectrum measured from 
SDSS DR7 LRG sample has been released recently [i^, 
with 45 band power measurements at comoving scales 
k = 0.0221-0.199 /iMpc~^, corresponding band-power 
window functions, and inverse covariance matrix between 
measurement errors. For each GCG model of a-wxo, 
we obtain 45 data points by convolving the linear bary- 
onic matter power spectrum with the band-power win- 
dow functions and normalize the point at the largest 



7 



scale {k = 0.0221 /iMpc"^) to that of (convolved) fidu- 
cial ACDM model power spectrum based on the WMAP 
5-year observation. Then we estimate the GCG model 
probability distribution £ oc ' on a-wxo plane by 
scanning = d"^C~^d, where d is a 45 x 1 vector con- 
taining GCG powers relative to LRG measurement and 
C is the 45 x 45 covariance matrix. During the scan- 
ning, other cosmological parameters have been fixed. By 
adding extremely large noise on the diagonal components 
of covariance matrix corresponding to k > 0.1 /iMpc~^ 



scales, we have effectively excluded LRG power spectrum 
information at scales where non-linear clustering dom- 
inates. Notice that during the analysis we have com- 
pletely ignored non-linear clustering properties of the 
matter power spectrum, and thus have not derived the 
halo power spectrum to compare with measured LRG 
power spectrum. For the case of ACDM-motivated mock 
power spectrum data (Fig. [S]), similar analysis has been 
done but without convolution operations. 
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